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Abstract 

We study the regularity of the extremal solution of the semilinear biharmonic equation A 2 u = rrjjjp , 
which models a simple Micro- Electromechanical System (MEMS) device on a ball B C ffi , under 
Dirichlet boundary conditions u = d v u — on dB. We complete here the results of F.H. Lin and Y.S. 
Yang [14] regarding the identification of a "pull-in voltage" A* > such that a stable classical solution 
it A with < u\ < 1 exists for AG (0, A*), while there is none of any kind when A > A*. Our main result 
asserts that the extremal solution u\* is regular (sup s u\* < 1) provided N < 8 while u\* is singular 
(sup s u\* — 1) for N > 17, in which case 1 — Co|a;| 4/ ' 3 < u\*(x) < 1 — 2?| 4//3 on the unit ball, where 

C := I f j' and A := The singular character of the extremal solution for the remaining 

cases (i.e., when 9 < N < 16) requires a computer assisted proof and will not be addressed in this paper. 



1 Introduction 

The following model has been proposed for the description of the steady-state of a simple Electrostatic 
MEMS device: 

aA 2 u = ((3 f \Vu\ 2 dx + 7) Au + — r in ft 

< u < 1 in O W 

u = ad u ii = on dQ, 

where a, /3,7, X ^ 0, / e C(S7, [0, 1]) are fixed, Q is a bounded domain in R N and A > is a varying 
parameter (see for example Bernstein and Pelesko 19 ). The function u(x) denotes the height above a point 
x G 57 C M. N of a dielectric membrane clamped on dCl, once it deflects forwards a ground plate fixed at 
height z = l, whenever a positive voltage - proportional to A - is applied. 

In studying this problem, one typically makes various simplifying assumptions on the parameters a, /3, 7, x, 
and the first approximation of ([T]) that has been studied extensively so far is the equation 

-Au = XJHj in n 

< u < 1 in fl (S)x,f 

u = on d£l, 
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where we have set a = [3 = \ = and 7=1 (see for example [51 [HI IS] and the monograph 7\) . This simple 
model, which lends itself to the vast literature on second order semilinear eigenvalue problems, is already 
a rich source of interesting mathematical problems. The case when the "permittivity profile" / is constant 
(/ = 1) on a general domain was studied in [IB], following the pioneering work of Joseph and Lundgren 
[15] who had considered the radially symmetric case. The case for a non constant permittivity profile / was 
advocated by Pelesko [TH] , taken up by [TT] , and studied in depth in [BJ [HI [3] • The starting point of the 
analysis is the existence of a pull-in voltage A*(f2, /), defined as 



It is then shown that for every < A < A*, there exists a smooth minimal (smallest) solution of (S)\j, while 
for A > A* there is no solution even in a weak sense. Moreover, the branch A h- > u\(x) is increasing for each 
x £ f2, and therefore the function u*(x) := lim\y\* u\(x) can be considered as a generalized solution that 
corresponds to the pull-in voltage A*. Now the issue of the regularity of this extremal solution - which, by 
elliptic regularity theory, is equivalent to whether sup n u* < 1 - is an important question for many reasons, 
not the least of which being the fact that it decides whether the set of solutions stops there, or whether a new 
branch of solutions emanates from a bifurcation state (u*,A*). This issue turned out to depend closely on 
the dimension and on the permittivity profile /. Indeed, it was shown in [9] that u* is regular in dimensions 
1 < N < 7, while it is not necessarily the case for N > 8. In other words, the dimension N = 7 is critical 
for equation (S)\ (when / = 1, we simplify the notation (5')a,i into (S)\). On the other hand, it is shown 
in [5] that the regularity of u* can be restored in any dimension, provided we allow for a power law profile 
\x\ v with 7/ large enough. 

The case where (3 = 7 = \ = (and a = 1) in the above model, that is when we are dealing with the 
following fourth order analog of (S)\ 



was also considered by [H [14] but with limited success. One of the reasons is the lack of a "maximum 
principle" which plays such a crucial role in developing the theory for the Laplacian. Indeed, it is a well 
known fact that such a principle does not normally hold for general domains f2 (at least for the clamped 
boundary conditions u = d v u = on dft) unless one restricts attention to the unit ball = B in t*, 
where one can exploit a positivity preserving property of A 2 due to T. Boggio [3]. This is precisely what 
was done in the references mentioned above, where a theory of the minimal branch associated with (P)\ is 
developed along the same lines as for (S)\. The second obstacle is the well-known difficulty of extracting 
energy estimates for solutions of fourth order problems from their stability properties. This means that the 
methods used to analyze the regularity of the extremal solution for (S)\ could not carry to the corresponding 
problem for (P)\. 

This is the question we address in this paper as we eventually show the following result. 

Theorem 1.1. The unique extremal solution u* for (-P)a* in B is regular in dimension 1 < N < 8, while it 
is singular (i.e, sup B u* — 1) for N > 17. 

Actually, we believe that the critical dimension for (P)\ in B is TV = 8, as opposed to being equal to 7 in 
(S)\. We add that our methods are heavily inspired by the recent paper of Davila et al. [5] where it is shown 
that N = 12 is the critical dimension for the fourth order nonlinear eigenvalue problem 



A*(f2, /) := sup < A > : there exists a classical solution of (S)\j 




(P)x 




in B 
on dB, 



while the critical dimension for its second order counterpart (i.e., the Gelfand problem) is iV = 9. 
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Throughout this paper, we will always consider problem (P)\ on the unit ball B. We start by recalling some 
of the results from [4] concerning (P)\, that will be needed in the sequel. We define 

A* := sup | A > : there exists a classical solution of (P)a|, 

and note that we are not restricting our attention to radial solutions. We will deal also with weak solutions: 
Definition 1.1. We say that u is a weak solution of (P)\ if < u < 1 a. e. in B, ^\yi £ L 1 {B) and 



f uA 2 <t> = X [ — — — , G C\B) n H 2 {B) 

J B JB U _ u ) 



We say that u is a weak super- solution (resp. weak sub-solution) of (P)\ if the equality is replaced with the 
inequality > (resp. <) for all <j) G C 4 {B) n H$(B) with <f> > 0. 

We also introduce notions of regularity and stability. 

Definition 1.2. Say that a weak solution u of {P)\ is regular (resp. singular) if ||u||oo < 1 (resp. =) and 
stable (resp. semi-stable) if 

firiu) = inf {^(A^) 2 - 2X J b J^^ : <t> e Hg(B), \\<f>\\v = 1 
is positive (resp. non-negative). 

The following extension of Boggio's principle will be frequently used in the sequel (see Lemma 16] and 
[SJ Lemma 2.4]): 

Lemma 1.1 (Boggio's Principle). Let u G i 1 (B). Then u > a.e. in B, provided one of the following 
conditions hold: 

1. tie C 4 (B), A 2 u > on B, and u = = on dB. 

2. J B uA 2 (/)dx > for all < <j) G C 4 (B) n H$(B). 

3. ue H 2 (B), u = and < on dB, and J B AuA0 > for all < <j) G H 2 (B). 
Moreover, either u = or u > a.e. in B. 

The following theorem summarizes the main results in [1] that will be needed in the sequel: 
Theorem 1.2. The following assertions hold: 

1. For each < A < A* there exists a classical minimal solution u\ of (P)\. Moreover u\ is radial and 
radially decreasing. 

2. For A > A*, there are no weak solutions of (P)\. 

3. For each x G B the map A t— > u\(x) is strictly increasing on (0, A*). 

4- The pull-in voltage A* satisfies the following bounds: 

f32(10JV- N 2 - 12) 128 - 240iV + 72N 2 1 _ 4i/ x 

max <^ — -, } < A < 

\ 27 ' 81 J " " 27 

where v\ denotes the first eigenvalue of A 2 in H 2 (B). 

5. For each < A < A* 7 u\ is a stable solution (i.e., Hi(u\) > 0). 
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Using the stability of u\, it can be shown that u\ is uniformly bounded in Hq(B) and that 1 2 Ux is uniformly 
bounded in L 3 (B). Since now A i— > u\(x) is increasing, the function u*(x) := Mxo.x/'x* u\(x) is well defined 
(in the pointwise sense), u* G H§(B), 1 _ 1 M , G L 3 (B) and w* is a weak solution of (P)a« ■ Moreover u* is the 
unique weak solution of (-P)a* • 

The second result we list from [4] is critical in identifying the extremal solution. 

Theorem 1.3. If u G Hq(B) is a singular weak solution of (P)\, then u is semi-stable if and only if 
(u,A) = (u*,A*). 

2 The effect of boundary conditions on the pull-in voltage 

As in [5], we are led to examine problem {P)\ with non-homogeneous boundary conditions such as 

A 2 u = n A v> in fl 

a< u < 1 in B (P)x, a ,p 

u = a , c^it = (3 on 3_B, 

where a, [3 are given. 

Notice first that some restrictions on a and (3 are necessary. Indeed, letting &(x) := (a — f ) + f |a;| 2 denote 
the unique solution of 

A 2 $ = in B . , 

$ = a , = /3 on 95, ^ ^ 

we infer immediately from Lemma 1 1.1 1 that the function u — $ is positive in B, which yields to 

sup $ < sup u < 1 . 

B B 

To insure that $ is a classical sub-solution of (P)\, a ,p, w e impose a / 1 and (3 < 0, and condition sup$ < 1 
rewrites as a - | < 1. We will then say that the pair (a, f3) is admissible if /3 < 0, and a — ^ < 1. 

This section will be devoted to obtaining results for {P)\, a ,p when (a,/?) is an admissible pair, which are 
analogous to those for (P)\. To cut down on notation, we shall sometimes drop a and (3 from our expressions 
whenever such an emphasis is not needed. For example in this section u\ and u* will denote the minimal 
and extremal solution of (P)\ ta ,p- 

We now introduce a notion of weak solution for (P)\, a ,p- 

Definition 2.1. We say that u is a weak solution of (P)\ t a,p if a < u < 1 a.e. in B, rpz^p S ^(B) and if 

f(u- $)A 2 = A / 7 -^-, G C\B) n # 2 (B), 

where $ is given in (^). We say that u is a weak super- solution (resp. weak sub-solution) of (P)\, a ,p if the 
equality is replaced with the inequality > (resp. <) for f> > 0. 

We now define as before 

A* := sup{A > : (P)\, a ,p has a classical solution} 

and 

A* := sup{A > : (P)\. a ,p has a weak solution}. 

Observe that by the Implicit Function Theorem, one can always solve (P)\, a ,p for small A's. Therefore, A* 
(and also A*) is well defined. 

Let now U be a weak super-solution of (P)\, a ,i3- Recall the following standard existence result. 
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Theorem 2.1 (|2J). For every < / G L ' (B), there exists a unique < u G ^(B) which satisfies 

uA 2 c/ ) = f ft/, 
Jb Jb 

for all4>&C A (B)C\H%{B). 

We can now introduce the following "weak iterative scheme": Start with uq = U and (inductively) let u n , 
n > 1, be the solution of 

/ K-$)A 2 = A / v^GC 4 (B)ni/ 2 (i?) 

Jb (1 - -"„-ir 

given by Thcorcm l2.ll Since is a sub-solution of (P)x,a,B, one can easily show inductively by using Lemma 
11.11 that a < it n +i < u n < U for every n > 0. Since 

(i- u „)- 2 < (i-U)- 2 eL\B), 

we get by Lebesgue Theorem, that the function u = lim u n is a weak solution of (P)\ a g such that 

n — >+oo 

a < u < U . In other words, the following result holds. 

Theorem 2.2. Assume the existence of a weak super- solution U of (P)\ t a,p- Then there exists a weak 
solution u of (P)x,a,B so that a < u < U a.e. in B. 

In particular, we can find a weak solution of (P)\, a ,6 for every A G (0, A*). Now we show that this is still 
true for regular weak solutions. 

Theorem 2.3. Let (a,0) be an admissible pair and let u be a weak solution of (P)\, a ,f3- Then for every 
< /i < A, there is a regular solution for (P) ^a^p ■ 

Proof. Let e £ (0, 1) be given and let u = (1 — e)u + e$, where $ is given in (J2j) . We have that 

supu < (1 — e) + esup$ < 1 , inf u > (1 — e)a + einf $ = a, 
b b B B 

and for every < <j) G C 4 (B) n HqC 8 ) there holds: 



(1 -s + e($- l)) 2 " v ' J B (l-u) 2 ' 

Note that < (1 — s)(l — u) = 1 — u + e($ — 1) < 1 — u. So u is a weak super-solution of {P)(\- e ) 3 \,a,p 
satisfying supw < 1. From Theorem 12.21 we get the existence of a weak solution w of (P)h— E )3\ a a so that 

B 

a < w < u. In particular, supu> < 1 and w is a regular weak solution. Since s G (0, 1) is arbitrarily chosen, 

B 

the proof is complete. □ 



Theorem 12 .31 implies in particular the existence of a regular weak solution U\ for every A G (0, A*). Introduce 



now a "classical" iterative scheme: uq = and (inductively) u n = v n + n > 1, where v„ G H 2 (B) is the 



(radial) solution of 

A 2 «„ = A 2 ( Un -$) = - in B. (3) 

(1 - U„-l) 2 

Since u„ G Hq(B), u n is also a weak solution of ([3]), and by Lemma \l. II we know that a < u n < it„+i < U\ 
for every n > 0. Since supu„ < sup [/a < 1 for n > 0, we get that (1 — u„_i)~ 2 G L 2 (B) and the 
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existence of v n is guaranteed. Since v n is easily seen to be uniformly bounded in Hq(B), we have that 
u\ :— lim u n does hold pointwise and weakly in H 2 (B). By Lebesgue Theorem, we have that u\ is a 

/( — >+oo 

radial weak solution of (P)\, a ,/3 so that sup«A < swpU\ < 1- By elliptic regularity theory [T] u\ G C°°(B) 

B B 

and u\ — $ = d v (u\ — $) = on dB. So we can integrate by parts to get 

/ AW= / A 2 (u x -<5>)4> = / (u x - $)A 2 = A / 
Jb Jb Jb (1 - «*r 

for every </> € C A (B) n Hq(B). Hence, u A is a radial classical solution of (P)a,<*,/3 showing that A* = A*. 
Moreover, since i> and v\ := u\ — $ arc radially decreasing in view of [20], we get that ma is radially 
decreasing too. Since the argument above shows that u\ <U for any other classical solution U of (P)^, a .j3 
with fx > A, we have that ma is exactly the minimal solution and u\ is strictly increasing as A j A*. In 
particular, we can define u* in the usual way: u*(x) = lim u\(x). 

Finally, we show the finiteness of the pull-in voltage. 

Theorem 2.4. If (a, (3) is an admissible pair, then X*(a,(3) < +oo. 

Proof. Let u be a classical solution of (P)\, a ,p and let (ip, v\) denote the first eigenpair of A 2 in H${B) with 
ip > 0. Now, let C be such that 

/ (/3 Aip - ad v Aip) =C [ ip. 
JdB Jb 

Multiplying (P)\, a ,p by ip an d then integrating by parts one arrives at 

/ fl ((T^-"- c )* = °- 

Since ip > there must exist a point x £ B where ( 1 _^ S ))2 — Viu(x) — C < 0. Since a < < 1, one can 
conclude that A < sup a<u<1 (i/]W + C)(l — u) 2 , which shows that A* < +oo. □ 

The following summarizes what we have shown so far. 

Theorem 2.5. If (a, (3) is an admissible pair, then A* G (0, +oo) and the following hold: 

1. For each < A < A* there exists a classical, minimal solution u\ of (P)x,a,0- Moreover u\ is radial 
and radially decreasing. 

2. For each x G B the map A h- > u\(x) is strictly increasing on (0, A*). 

3. For A > A* there are no weak solutions of (P)x,a,p- 



2.1 Stability of the minimal branch of solutions 

This section is devoted to the proof of the following stability result for minimal solutions. We shall need 
yet another notion of H 2 (B)— weak solutions, which is an intermediate class between classical and weak 
solutions. 

Definition 2.2. We say that u is a H 2 (B) — weak solution of (P)x,a,i3 «/u — $ G H 2 (B), a < u < 1 a.e. in 
B > TT^p eL^B) and if 

f AuAcj, = A / V</> G C\B) n H 2 (B), 

Jb Jb U — u ) 

where $ is given in {![). FFe say that u is a H 2 (B) — weak super- solution (resp. H 2 (B) — weak sub-solution) 
°f (P)\,a,/3 if for <j) > the equality is replaced with > (resp. <) and u > a (resp. <), d v u < (3 (resp. >) 
on dB. 
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Theorem 2.6. Suppose (a, (3) is an admissible pair. 

1. The minimal solution u\ is then stable and is the unique semi-stable H 2 (B) — weak solution of (P)x,a./3- 

2. The function u* := lim u\ is a well-defined semi-stable H 2 (B) — weak solution of (P)\*, a ,@- 

3. Whenu* is classical solution, then /j,\(u*) = and u* is the unique H 2 {B)~ weak solution of {P)\* iQj| g. 
4- If v is a singular, semi-stable H 2 (B) — weak solution of (P)\ <a ,0i then v — u* and A = A* 

The crucial tool is a comparison result which is valid exactly in this class of solutions. 

Lemma 2.1. Let (a,/3) be an admissible pair and u be a semi-stable H 2 (B) — weak solution of {P)\ >a ,p. 
Assume U is a H 2 (B) — weak super- solution of (P)\. a ,f3 so that U — $ G Hq(B). Then 

1. u < U a.e. in B; 

2. If u is a classical solution and Hi(u) = then U = u. 

Proof, (i) Define w := u—U. Then by the Moreau decomposition [T7] for the biharmonic operator, there exist 
w\ 1 W2 G H 2 {B), with w = W\ + u> 2 , W\ > a.e., A 2 u>2 < in the H 2 (B)— weak sense and J B Au>iAu> 2 = 0. 
By Lemma 11.1 [ we have that W2 < a.e. in B. 
Given now < <f> € C£°(B), we have that 

/ AwAtj) < A / (/(«) - f(U))<f>, 

JB J B 

where f(u) := (1 — u)~ 2 . Since u is semi-stable, one has 

A / f\u)w\ < f (A Wl ) 2 = f AwA Wl < A / (f(u)-f(U)) Wl . 



Since w± > w one also has 



which once re-arranged gives 



Jj'{u)w Wl < j (flu) - /{Our,. 



fW! > 0, 



where f(u) — f(u) — f(U) — f'(u)(u — U). The strict convexity of / gives / < and / < whenever u ^ U. 
Since Wi > a.e. in B one sees that w < a.e. in B. The inequality u < U a.e. in B is then established. 



ii) Since u is a classical solution, it is easy to see that the infimum in is attained at some (j). The 

2A 
(l-u) 3 



function <f) is then the first eigenfunction of A 2 — 2A , 3 in Hq(B). Now we show that <f> is of fixed sign. 



Using the above decomposition, one has <f> — (j>\ + 02 where <pi G Hq(B) for i = 1, 2, 0i > 0, f B AfaAfo = 
and A 2 02 < in the Hq(B)— weak sense. If <j> changes sign, then <j)\ ^ and 02 < in B (recall that either 
4>2 < or 02 = a.e. in B). We can write now: 

/ B (A(0! - 2 )) 2 - Af (uK0i - 2 ) 2 ^ / B (A0) 2 - Af( M )0 2 

in view of 0102 < — 0102 in a set of positive measure, leading to a contradiction. 

So we can assume > 0, and by the Boggi's principle we have > in B. For < t < 1 define 

g{t)= j A [tU + (1 - t)u] A0 - A / f(tU + (l-t)u)(f), 
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where <b is the above first eigenfunction. Since / is convex one sees that 

g{t) > A / [tf(U) + (l - t)f(u) - f(tu + (l - t)u)] 4> > o 
Jb 

for every t > 0. Since 5(0) = and 

</(0) = / A(U-u)A<t>-\f'(u)(U-u)<t> = 0, 



we get that 

g"(0) = -\ I f'(u)(u -u) 2 <b>Q. 



Since f"(u)4> > in i?, we finally get that U = u a.e. in £?. □ 
Based again on Lemma [Ol3), we can show a more general version of the above Lemma [2TT1 

Lemma 2.2. Let (a,/3) be an admissible pair and /3' < 0. Let u be a semi-stable H 2 (B) — weak sub-solution 
°f (P)x a 8 with u = a, d y u — (3' > (3 on dB. Assume that U is a H 2 (B) — weak super- solution of (P)\,a 8 
with U = a, d v U — P on dB. Then U > u a.e. in B. 

Proof. Let u 6 H 2 (B) denote a weak solution to A 2 u — A 2 (it — U) in B. Since u — u + U = and 
d v (u — u + U) < on dB, by Lemma 11.11 one has that u > u — U a.e. in B. Again by the Moreau 
decomposition [17] , we may write u as u = w + v, where w,v G Hq(B), w > a.e. in B, A 2 v < in a 
H 2 (B)-weak sense and J B AwAv = 0. Then for < 4> £ C 4 (B) n H%(B) one has 



AuA(/)= / A(u - U)A(b < A / (f(u)-f(U))(f>. 
Jb Jb 

f AiiAw <\f (f(u) - f(U))w. 
Jb Jb 

A f f'{u)w 2 < [ {Aw) 2 = [ AiAw, 
Jb Jb Jb 



In particular, we have that 
Since by semi-stability of u 

we get that 

/ f'(u)w 2 < f (f(u)-f(U))w. 

J B J B 

By Lemma ll.ll we have v < and then w > u > u — U a.e. in B. So we see that 

0< / (/(«)- /(CO -/'(«)(« -U))w. 



The strict convexity of / implies as in Lemma I2TT1 that U > u a.e. in B. □ 
We shall need the following a-priori estimates along the minimal branch u\. 
Lemma 2.3. Let (a, j3) be an admissible pair. Then one has 

(u x - $) 2 f u x - $ 



(l-O 3 - Jb (1-«a) 2 ' 
where $ is given in 0). In particular, there is a constant C > so that for every A 6 (0, A*), we have 
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Proof. Testing (P)\, a ,p on u x - $ G C 4 {B) n flj 2 (5), we sec that 



A / 7^^ = / A« A A(« A -$)= / (A(u A -$)) 2 >2A / 
7b (1 - u A ) 2 7 b Jb (1 



(1 - u x ) 2 J B ' 7b 7b (1 ~ u x ) 3 

in view of A 2 <E> = 0. In particular, for 6 > small we have that 

1 < 1 / K - <&) 2 1 /■ 1 



{\u x -$\>8} (1 - "a) 3 <5 2 7{|„ A -*|>5} (1 - w A ) 3 (5 2 7b (1 - "a) 2 

7{|n A -#|>5} (1 - Ux) 3 

by means of Young's inequality. Since for <5 small, 

1 



(|«a-*|<<5} i 1 _ ""a) 3 

for some C > 0, we can deduce that for every A G (0, A*), 

1 



for some C > 0. By Young's and Holder's inequalities, we now have 

/ (Au x ) 2 = f Au x A<S> + A / " A ~ * < J / (Au A ) 2 + C s + c( f 
Jb Jb Jb I 1 - u \) Jb \Jb 



and estimate Q is therefore established. □ 



We are now ready to establish Theorem 12.61 

Proof (of Theorem I2.6[) : (1) Since |ua||oo < 1, the infimum defining Hi(u x ) is achieved at a first 
eigenfunction for every A G (0, A*). Since A i— » ua(:e) is increasing for every x G B, it is easily seen that 
A i ► fJ-i(u x ) is an increasing, continuous function on (0, A*). Define 

A** := sup{0 < A < A* : Hi{u x ) > 0}. 

We have that A** = A*. Indeed, otherwise we would have that fj,i(u x „) — 0, and for every n G (A**, A*) 
would be a classical super-solution of (P)a»,,q,/3- A contradiction arises since Lemma 12 . 1 1 implies — u x „- 
Finally, Lemma |2~T1 guarantees uniqueness in the class of semi-stable H 2 (B) — weak solutions. 

(2) By estimate (fj| it follows that u x — > u* in a pointwise sense and weakly in H 2 (B), and ^ . G L 3 (B). 
In particular, it* is a H 2 (B)— weak solution of (P)a*,q./3 which is also semi-stable as limiting function of the 
semi-stable solutions {u x }. 

(3) Whenever ||m*||oo < 1, the function u* is a classical solution, and by the Implicit Function Theorem we 
have that = to prevent the continuation of the minimal branch beyond A*. By Lemma 12.11 u* is 
then the unique H 2 (B)— weak solution of (P)\* t a,/3- An alternative approach -which we do not pursue here- 
based on the very definition of the extremal solution u* is available in [4] when a = (3 = (see also [15]) to 
show that u* is the unique weak solution of (P) x *, regardless of whether u* is regular or not. 

(4) If A < A*, by uniqueness v = u x . So v is not singular and a contradiction arises. 

By Theorem l2.5f 3') we have that A = A*. Since v is a semi-stable H 2 (B)— weak solution of (P) x * >a ,i3 and u* 
is a H 2 (B)— weak super-solution of {P)\* t a,p-, we can apply Lemma [2~T1 to get v < u* a.e. in B. Since u* is 
a semi-stable solution too, we can reverse the roles of v and u* in Lemma |2. II to see that v > u* a.e. in B. 
So equality v = u* holds and the proof is done. 
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3 Regularity of the extremal solution for 1 < N < 8 

We now return to the issue of the regularity of the extremal solution in problem (P)\. Unless stated 
otherwise, u\ and u* refer to the minimal and extremal solutions of (P)\- We shall show that the extremal 
solution u* is regular provided 1 < N < 8. We first begin by showing that it is indeed the case in small 
dimensions: 

Theorem 3.1. u* is regular in dimensions 1 < N < 4. 

Proof. As already observed, estimate Q implies that f(u*) = (1 — u*)~ 2 G L^(B). Since u* is radial and 
radially decreasing, we need to show that u*(0) < 1 to get the regularity of u* . The intcgrability of f(u*) 
along with elliptic regularity theory shows that u* G W 4 '? (B). By the Sobolev imbedding Theorem we get 
that u* is a Lipschitz function in B. 
Now suppose u*(0) = 1 and 1 < N < 3. Since 

1 C 

> T-r in B 



1 — u \x\ 

for some C > 0, one sees that 



= f' I i < / ^-^3 < 00. 



B 



x\ 3 ~ J B (l-O 



A contradiction arises and hence u* is regular for 1 < N < 3. 

For N — 4 we need to be more careful and observe that u* S C x '3(-B) by the Sobolev Imbedding Theorem. 
If u*(0) = 1, then Vw*(0) = and 

1 C . „ 
> — — in B 

for some C > 0. We now obtain a contradiction exactly as above. □ 
We now tackle the regularity of u* for 5 < N < 8. We start with the following crucial result: 
Theorem 3.2. Let N > 5 and (u* , A*) 6e t/ie extremal pair of (P)\- When u* is singular, then 

1 -u*(x) < C \x\i inB, 

where C := (f )* and A := gi^^i) , 

Proof. First note that Theorem 1 1.2( 4) gives the lower bound: 

i28~ 240iV + 72iV 2 
^ A = si • (5) 

i 

For <5 > 0, we define ug(x) := 1 — C,j|a;|3 with Cs '■= (x^^) 3 > Since N > 5, we have that 
u & e ff( 2 „ c ( RiV ). TXI7 e ^L^) and u « is a -ff 2 -weak solution of 

A n A* + 5A w 

(1 - "i) 2 

We claim that 14,5 < u* in B, which will finish the proof by just letting 8 — ► 0. 

Assume by contradiction that the set T := {r G (0,1) : u$(r) > u*(r)} is non-empty, and let r\ = sup T. 
Since 

u i (l) = l-Ca<0 = u*(l), 
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we have that < r\ < 1 and one infers that 

a := u*(ri) = Us(n) , := (u*)'(ri) > u' s (ri). 

_ 4 

Setting us t n (f) = r x 3 (us(r±r) — 1)+1, we easily see that us,n is a H 2 (B) — weak super-solution of (P)x*+5Xa' B'j 
where 

a' : = r ^ 3 (a - 1) + 1 , [3':=r^(3. 

_ 4 

Similarly, let us define it* (r) = r 1 3 (u*(rir) — 1) + 1. The dilation map 

_ 4 

w — > w ri (r) = r x 3 (iy(nr) — 1) + 1 (6) 
is a correspondence between solutions of (P)a on P and of (P) 4 on B -i which preserves the 

P 2 — integrability. In particular, (it* , A*) is the extremal pair of (P) _4 on B 1 -i (defined in the 

X,l— r t 3 ,0 1 

obvious way). Moreover, u* is a singular semi-stable H 2 {B)— weak solution of (P) x* ,a' ,/3' ■ 

Since it* is radially decreasing, we have that /?' < 0. Define the function w as ui(x) := (a' — ^-) + ^-|x| 2 +7(x), 
where 7 is a solution of A 2 7 = A* in B with 7 = <9„7 = on dB. Then w is a classical solution of 

A 2 w = A* in B 

w = a' , d v w — (3' on dB. 



Since p > A*, by Lemma [TTT1 we have u* ± > w a.e. in B. Since w(0) = a' — + 7(0) and 7(0) > 0, 

the bound u* x < 1 a.e. in B yields to a' — 4- < 1. Namely, (a', ft') is an admissible pair and by Theorem 
2.6f 4) we get that (u* l5 A*) coincides with the extremal pair of (P)x,a',3' in P- 



Since (a',/3') is an admissible pair and is a H 2 (B)— weak super-solution of (P)a*+5a q ' a', by Theorem 
12 .21 we get the existence of a weak solution of (P)a*+5A a' fl'- Since A* + 5X > A*, we contradict the fact that 
A* is the extremal parameter of (P)x,a',B'- d 

Thanks to this lower estimate on it*, we get the following result. 

Theorem 3.3. If5<N<8, then the extremal solution u* of (P)a is regular. 

Proof. Assume that u* is singular. For e > set ip(x) := \x\~ he and note that 

N 2 ( N - A) 2 

(AVO 2 = (H N + 0(e))\x\- N + 2 ^ where H N := 1 - 1 . 

16 

Given ij £ Cq°(B), and since iV > 5, we can use the test function nip £ H 2 (B) into the stability inequality 
to obtain 

2A / n ^3 g / (A^) 2 + 0(i), 

where 0(1) is a bounded function as e \ 0. By Theorem 13.21 we find that 
and then 

2A f \x\~ N+2e < (H N + 0(s)) f \x\- n+2e + 0(1). 

JB J B 

Computing the integrals one arrives at 

2A < H N + O(e). 

As e — > finally we obtain 2A < H^. Graphing this relation one sees that N > 9. □ 
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We can now slightly improve the lower bound ((SJ). 



Corollary 3.1. In any dimension N > 1, we have 

A > A = 3ii ( 7 ) 

9 

Proof. The function u := 1 — |a;|$ is a H 2 (B)— weak solution of (P)a 0-4- ^ ^y contradiction A* = A, then 
u is a H 2 {B)— weak super-solution of (P)a for every A € (0, A*). By Lemma |2. II we get that u\ < u for all 
A < A*, and then u* < u a.e. in B. 

If 1 < N < 8, u* is then regular by Theorems 1331 and [331 By Theorem [2^ 3) there holds fJ,i(u*) = 0. 
Lemma 12.11 then yields that u* — u, which is a contradiction since then u* will not satisfy the boundary 
conditions. 



If now N > 9 and A = A*, then Co = 1 in Theorem 13.21 and we then have u* > u. It means again that 
u* = u, a contradiction that completes the proof. □ 



4 The extremal solution is singular for N > 17 

In this section, we will need the following improved Hardy-Rellich inequality, which is valid for N > 5 (see 
[10] and references therein): 

/ (A^) 2 >H N [ f^ + C [ V^e Hl(B) 

Jb Jb \ x \ Jb 

where Pjv := — — — is optimal and C > 0. As in the previous section (u* , A*) denotes the extremal pair 
of (P)a- We first show the following upper bound on u* . 

Lemma 4.1. If N >9, then u* < 1 — \x\i in B. 

Proof. Recall that A := ° v 3 ^ 3> < A*. If A = A*, then by the proof of Corollary O we know that 
u* < u. 

Suppose now that A < A*. We claim that u\ < u for all A € (A, A*). Indeed, fix A and assume by contradiction 
that 

Pi := inf{0 < P < 1 : u x < u in (P, 1)} > 0. 

From the boundary conditions, one has that u\{r) < u(r) as r — > 1~. Hence, < Pi < 1, a := u\(R\) — 
u(Ri) and /3 :— w' A (Pi) < u'(Ri). Introduce, as in the proof of Theorem [321 the functions (ux)^ and 
We have that (ux)^ is a classical super-solution of (P)A,a',/3'j where 

a' := R^ {a - 1) + 1 , 0' := [3. 

Note that (u)_R d is a H 2 (B)— weak sub-solution of (P)a a ' B' which is also semi-stable in view of 2A < Pjv 
and the Hardy inequality. By Lemma |2.2[ we deduce that {ux)^ > (u)b.i i n B. Note that, arguing as in 
the proof of Theorem 13. 2[ (a',/3') is an admissible pair. 

We have therefore shown that ux > u in B^ and a contradiction arises in view of the fact that lim u(x) = 1 
and ||ua||oo < 1- It follows that ux < u in B for every A 6 (A, A*), and in particular u* < u in B. □ 

Our approach for showing that u* is singular for large dimensions, will depend on the sign of Hn — 2A*. 

Theorem 4.1. If N >9 and A* < f/ien i/ie extremal solution u* of (P)x is singular. 
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Proof. Let ip € C^°(B) with J B ip 2 — 1. By Lemma FPl and the improved Hardy-Rellich inequality (see jlOj). 
one then has 

It follows that > and u* must be singular, since otherwise, one could use the Implicit Function 

Theorem to continue the minimal branch beyond A* . □ 

We can now show the following result about the extremal solution. 
Theorem 4.2. The following upper bounds on A* hold in large dimensions. 

1. IfN > 31, then A* < 27A < 4f . 

2. If 17 < N < 30, then A* < 

The extremal solution is therefore singular for dimension N > 17. 

Proof. Consider for any m > the following function: 

w m := 1 - 3m/(3m - 4)r 4/3 + 4r m /(3m - 4) (8) 

Assume first that AT > 31, then 27A < We shall show that u>2 is a singular H 2 {B)— weak sub-solution 
of (P)27\ that is semi-stable. Indeed, write 

w 2 := 1 - \x\i - 2(\x\§ - \x\ 2 ) = u-<t> , 

where O := 2(|a;|i - |x| 2 ), and note that u> 2 £ Hq{B), G L 3 (B), < w 2 < 1 in S, and 

27A 

AV< 7 . ^ inB\{0}. 

(1 - W2T 

So W2 is H 2 (B)~ weak sub-solution of (-P) 2 7a- Moreover, since 27A < 4f , and since <^>o > 0, we get that 

54A / ^ — Bn f f <H N I n4 < / (AVO 2 

for all if> € Hq(B). Hence W2 is also semi-stable. If now 27A < A*, then by Lemma f2.2[ W2 is necessarily 
below the minimal solution u 2 -j\ which contradicts the fact that u>2 is singular. Hence A* < 27A < 

Now consider the function 

9 4 4 , 

11)3 '■= 1 r:i H — r ■ 

5 5 

We show that it is a singular H 2 {B)— weak sub-solution of (P) h n that is semi-stable. Indeed, we clearly 

2 

have that < w 3 < 1 a.e. in £?, W3 s Hq(B) and y^r/- G L 3 (B). To show the stability condition, we 
consider %j) S C^°(B) and write 

[ V [ib 2 I f ib 2 

H * / M A3 = 125i ^ / ~ < 125H N SUp j—— / X_ 

7b (1-103) 7s (9ra - 4r 3 ) 3 o<r<i (9 - 4r 3 ~3 )3 J B r 4 

JB r JB 

An easy computation shows that 



H N a2 25H n 9A 12 

A w 3 



2(1 -u> 3 ) 2 " 2(9rt-4r 3 )2 5 r f 5 r 

25N 2 {N-A) 2 8{N- |)(JV- j) 12 AT 2 - 1 
32(9rl - 4r 3 ) 2 5r§ 5 r 
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and by using Maple, one can verify that this final quantity is nonnegative on (0, 1), whenever 17 < N < 30, 
hence u>3 is a subsolution of (P) h n . If now, < A*, then Lemma [2.21 would imply that the minimal 

solution u\ is larger than and hence is singular for < A < A* , which is a contradiction. □ 

Remark 4.1. We believe that the extremal solution is singular for all N > 9 and for that, one needs to 
construct again for the remaining cases 9 < N < 16 (and ^f- < A*), a singular H 2 (B)— weak sub-solution 
of (P) h n that is semi-stable. However, one can show that (at least for N = 9) such a sub-solution cannot 

2 

be obtained by simply perturbing u with a function of the form cf>o = |/3r"(l — r' 3 ). 

The construction of such a sub-solution for the remaining cases, i.e. when 9 < N < 16 and < A*, will 
therefore very likely require a computer assisted proof, that we leave open to the interested reader. 
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